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ABSTRACT 



Of 



We present F850LP-F160W color gradients for 11 early-type galaxies (ETGs) at 
1.0<z spec <1.9 selected from the GOODS South field. Significant negative F850LP- 
F160W color gradients (core redder than the outskirts) have been detected in ~ 70% 
of our sample within the effective radius R e , the remaining 30% having a flat color 
profile consisten with a null gradient. Extending the analysis at R>R e , enclosing the 
whole galaxy, we have found that the fraction of high-z ETGs with negative F850LP- 
F160W color gradients rises up to 100%. For each galaxy, we investigate the origin 
of the radial color variation with an innovative technique based on the matching of 
both the spatially resolved color and the global spectral energy distribution (SED) to 
predictions of composite stellar population models. In fact, we find that the age of the 
stellar populations is the only parameter whose radial variation alone can fully account 
for the observed color gradients and global SEDs for half of the galaxies in our sample 
(6 ETGs), without the need of radial variation of any other stellar population prop- 
erty. For four out of these six ETGs, a pure metallicity variation can also reproduce 
the detected color gradients. Nonetheless, a minor contribution to the observed color 
gradients from radial variation of star-formation time scale, abundance of low-to-high 
mass stars and dust cannot be completely ruled out. For the remaining half of the 
sample, our analysis suggests a more complex scenario whereby more properties of the 
stellar populations need to simultaneously vary, likely with comparable weights, to 
generate the observed color gradients and global SED. Our results show that, despite 
the young mean age of our galaxies (<3-4 Gyr), they already exhibit significant differ- 
ences among their stellar content. We have discussed our results within the framework 
of the widest accepted scenarios of galaxy formation and conclude that none of them 
can satisfactorily account for the observed distribution of color gradients and for the 
spatially resolved content of high-z ETGs. Our results suggest that the distribution of 
color gradients may be due to ETGs forming by different mechanisms. 

Key words: galaxies: elliptical and lenticular, cD galaxies: evolution galaxies: 
formation galaxies: high-redshift galaxies: stellar content. 



1 INTRODUCTION 

A viable way to gather insight on the processes that concur 
to accrete the stellar mass in early-type galaxies (ellipticals 
plus SO's, hereafter ETGs) is to analyze the spatial distri- 
bution and properties of their stellar content which, in prin- 
ciple, can be directly connected to the events experienced 
by the galaxies. Indeed, the different scenarios proposed to 
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explain the formation of ETGs give different predictions on 
their stellar population content. 

The revised monolithic model predicts that, in a cold 
dark matter framework, massive ETGs assemble the bulk 
of their mass at z> 2-3 through the merger of small sub- 
structures moving in a common potential well. This ini- 
tial collapse might be regulated by cold gas streams from 
the cosmological surro undings, also known as cold accre- 
tion (|Dekel et alj|2009l ). the latters becoming less important, 
for high-mass galaxies, at lower redshift. The subsequent 
evolution should be mainly characterized by the aging of 
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their stellar populations, with small new episodes of star- 
formation at z<l, related, e.g., to the capture of satellites 
dKatdll99ll : iKawatallioOll : iKobavashillioO 4: lMerlin fc Chiosil 
During the gravitational dissipative collapse, the 
metal-enriched gas should naturally flows towards the cen- 
ter of the galaxy, leading to an ETG with stellar populations 
more metal-rich in the center than in the external regions 
(i.e. negative metallicity gradient). Moreover, because of the 
deeper potential well, the star-formation is expected to last 
longer in the central than in the outer regions. This would 
lead to null or mildly positive age gradients, with stellar 
populatio ns in the center ~ 10 % younger than those in the 
outskirts (|Kobavashill2004T ). Nonetheless the effect of metal- 
licity variation on color profile should be the dominant one, 
thus in this scenario ETGs are generally expected with neg- 
ative color gradients. 

The competing formation scheme is the hierarchical sce- 
nario. Following the hierarchical assembly of cosmic struc- 
tures, ETGs are supposed to for m through gas-rich ( " wet") 
mergers of disc ga laxies (e.g. iToomre fc Toomra 1 19721 ; 
iDe Lucia et~ail 2006) at high-redshift (z~4-5). In this phase, 
a large fraction of the stellar mass of a galaxy is assem- 
bled through central intense bursts of star formation (e.g. 
Rcnzini 2006). Concurrently, a "dry" merger picture has also 
been advocated, where bright ETGs would form thro ugh the 
merging of quiescent galaxies (e.g. iBell et alj|2004 ). Actu- 
ally, in the last years a new scheme of mass accretion of 
massive ETGs, known as inside-out growth, is becoming 
widely accepted. This scenario is motivated by the obser- 
vational evidence of ETGs at 1.0<z<2.5 with effective radii 
3-5 times smaller than t he mean radius of l ocal ETGs with 
the same stellar mass (e.g lDaddi et al.ll2005l ;l Longhetti et al.l 
120071 ). In thi s context, supported by a wealth o f sim- 



ulations (e.g. iKhochfar fc Sil 
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IWuvts et alTl2010l ; iNaab et al 
ETGs are supposed to be formed at high redshift (z~4-5) 
as compact spheroids result of gas-rich mergers. Then, at 
lower redshift, compact ETGs would undergo subsequent 
minor- "dry" mergers, whose main effect is to add an exter- 
nal low-mass density envelope to the compact core ETGs, 
enlarging the effective radius while leaving the stellar mass 
nearly constant. Indeed, this scenario shows some limita- 
tions, such as the not plausible number of minor- "dry" merg- 
ers necessary to enlarge the ETGs' size which could produce 
a scatter in the fundamental plane much larger than the ob- 
served one, and its failure to explain the presence of nor- 
mal ETGs obs erved at high-z in number similar to the com - 
pact ones (e.g. lNipoti et al.ll2009l : rSaracco et ai1l201ll , l2010h . 
From a theoretical point of view, the "wet"-merger scenario 
predicts ETGs with significant radial age variations. Indeed, 
the galaxy remnant of a wet merger should be characterized 
by a central stellar population younger, and more metal- 
rich, than the outer o ne, i.e. by a pos itive (negative) age 
(metallicity) gradient (|Kobav ashi 2004). On the contrary, 
"dry"-mergers, mixing the pre-existing stellar populations 
of progenitor galaxies, should dilute any radial variation of 
age and metallicity, produc ing flatter distributions (but see 
also lDi Matteo et all (|200gh ). In the inside-out scenario mi- 
nor dry-mergers are actually believed to only add an outer 
low-density envelope on top of a compact core, without mix- 
ing the pre-existing stellar-populations but redistributing 
the star content of the satellites in the outer regions of the 



compact ETG. 

The picture described so far shows how important can be 
to spatially resolve the properties of the underlying stellar 
populations of ETGs in order to pinpoint their formation 
scenario. The most viable way to carry out information on 
the radial variation of the stellar content of a galaxy is to 
investigate its radial color variation, being the color of a stel- 
lar population tightly dependent on its age, star-formation 
time scale, metallicity, and (eventual) presence of dust. 
ETGs at high redshift (hereafter with "high redshift" we 
mean z >1) represent the best-suited benchmark to investi- 
gate and constrain the mechanisms driving the mass assem- 
bly of spheroids due to the short time elapsed since their for- 
mation. So far, instrumental limits have allowed only studies 
of the global (i.e. integrated) properties of high-z spheroids, 
preventing, indeed, any measure of their spatially resolved 
information. Recently, the capabilities of the Hubble Space 
Telescope (HST) have made possible to overwhelm p art of 
these limitations for the first time. Gargiulo et al. (|201ll . 
hereafter GSL11) taking advantage of the deep and high- 
resolution HST Advanced Camera for Surveys (HST/ ACS) 
images of the GOODS South field derived F606W - F850LP 
color gradients (~ (UV -U) ,-estframe at z=1.5, A e // of F606W 
and F850LP filter ~ 5810 A and 9010 A respectively) for 20 
ETGs at 1 < z S pec < 2. In their work they ascertained the 
feasibility of this analysis up to z~2 and presented the first 
spatially resolved information for ETGs at high-z. Despite 
the short wavelength baseline covered, ~ 50 % of the galax- 
ies showed a significant (positive or negative) color gradi- 
ent. These results clearly showed that, after 3-4 Gyr (at 
most) from their birth, ETGs do not exhibit a unique spatial 
distribution of their stellar populations, implying that they 
followed different mass assembly paths. Unfortunately, the 
available HST data, mainly optical, prevented the authors to 
discriminate the drivers of the observed color gradients (e.g. 
radial variations of age, metallicity, ...) and consequently, to 
constrain the mechanisms responsible for them. Indeed, at 
<z> ~ 1.5, the galaxy emission sampled by the F606W and 
F850LP filters is sensitive to both age and dust variation. 
Moreover, it is dominated by the youngest (~ lGyr) stellar 
populations, missing any information on the distribution of 
the oldest stellar populations. 

The recent advent of the first HST Wide Field Camera 3 
(HST/WFC3) near infrared images for part of the GOODS 
South area (see Section 2 for details) has opened new possi- 
bilities to study color gradients of high-z ETGs and to con- 
strain their origin. In this paper we combine the information 
provided by the WFC3/F160W-band (~ R-band rest-frame 
at <z> ~ 1.5, Xeff of the filter ~ 15369 A) and F850LP- 
band emission to derive the F850LP-F160W color gradients 
(~ (U-R)restframe a t z=1.5) for a sample of 11 ETGs at 
<z> ~1.5. The bands we selected sample emissions domi- 
nated by different stellar populations. In particular, differ- 
ently from F606W-F850LP color, the F850LP-F160W color 
is much more sensitive to age variations than dust content. 
This allows us to extend and to complement the analysis 
presente d by GSL11 . 

Actually. IGuo et alj i|201ll) presented F850LP-F160W color 
gradients for 4 massive passively evolving galaxies in the 
GOODS South area at 1.3 < z 3pec < 2. They derived the 
color profiles by measuring optical-NIR colors in concentric 
annuli, and found that high-z ETGs have cores redder than 
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the outerskirts. The observed radial trend in color gradients 
is not reproduced by the radial variation of a single stel- 
lar population parameter (age, metallicity, extinction), al- 
though they found that dust should partially contribute to 
generate the observed color distribution. In this paper, we 
exami ne a sample three times larger than that of lGuo et al.l 
fioil]), deriving F850LP-F160W color profiles from the 2- 
dimensional fit of the light profiles in the two bands. Then, 
we present a new approach to exploit the wealth of available 
information to constrain the radial variation of the underly- 
ing stellar populations. In fact, taking advantage of a unique 
set of data and an innovative procedure, we are able to con- 
strain the radial variation of stellar population parameters 
(age, metallicity, dust, star-formation time scale, and initial 
mass function) and their contribution to produce the color 
gradients we observe in high-z ETGs. Finally, we compare 
our findings to the prediction of theoretical formation mod- 
els. 

Throughout the paper we adopt a standard ACDM cosmol- 
ogy with Ho=70kms _1 Mpc" 1 , f2 m =0.3 and fi A =0.7. All 
the magnitudes are in the AB system. 



2 THE SAMPLE 

We have derived F850LP-FI60W color gradients for If 
ETGs at <z> ~ f.5. The sample has been extracted from 
the complete sample of 34 ETGs selected on the whole 
GOO DS South field ( GOODS-South y2; G iavalisco et al. 
l2004h and presented in lSaracco et al.l l|2010l ). The morpho- 
logical classification was performed by the authors both on 
the basis of a visual inspection of the F850LP images and 
of Sersic index ufssolp {riFS50LP > 2). Starting from this 
complete sample we have restricted our preliminary analy- 
sis to those ETGs with near-IR WFC3 data (16 ETGs out 
34). Indeed, as pointed out in the Introduction, two areas 
of the GOODS South field have been imaged in the F160W 
band with the WFC3 extending the space-based data from 
the pre-existing optical domain to the near-infrared one. 
The Early Release Science (ERS, propID: 11359, PI: R. W. 
O'Connell) imaged an area of ~ 40 arcmin 2 . Briefly, the ob- 
servations were acquired with a 2-cycles-long exposures for 
a total exposure time of ~ 6 ks, reaching a 5<r d epth in the 
AB system of F160W = 27.25 (Windhorst et al. l201lh . Ad- 
ditionally, the HUDF09 HST Treasury program (GO 11563, 
PI: Illingworth) provides the first ultra-deep near-IR WFC3 
observation of the Hubble Ultra Deep Field (HUDF). The 
images cover an area of ~ 4.7 arcmin 2 . At the time of our 
analysis, the first epoch data were available for a total ex- 
posure time of ~ 80 ks and a depth at 5 a of 2 8.8 in the 
F160W filter (AB magnitude, Oesch et al. \201(j ). Starting 
from the raw WFC3 images, we have creat ed the mosaic im- 
ages with the software MULTIDRIZZLE dKoekemoer et all 
120021 ) reducing the pixel/scale of the mosaic from the orig- 
inal value of 0.128 "/px to 0.06 "/px. In the final mosaics 
we have stacked only those single exposures not affected 
by the presence of persistence. Thus, we have obtained the 
WFC3/F160W-band mosaics of the ERS and HUDF09 area. 
Both mosaics are characterized by a full width at half max- 
imum (FWHM) of ~ 0.2". The exposure times are ~ 80ks 
and ~ 6 ks for the HUDF09 and ERS, respectively. S/N > 5 
at 3r e are reached for the faintest galaxies of our sample in 



the ERS area, the one with lower exposure time. This S/N 
and the resolution of the WFC3 assure us reliable estimates 
of the surface brightness parameters even in the images with 
the shorter exposure time (see Sec. 3.1). 

Out of the 16 ETGs with NIR WFC3 data available, 
five objects have been excluded from the present analysis 
because of either potential problems in the modelling of the 
near-IR galaxy light distribution (4 ETG, Sec. 3), or of large 
uncertainties in the estimate of their structural parameters 
(1 ETG, Sec. 3.1). These selections result into a final sam- 
ple of 11 galaxies. Each one of these gala xies i s provided 
with spectroscopic redshift (Vanzella et al. 120081 and refer- 
ences therein) and with the wide photometric coverage of the 
GOODS survey: four deep HST/Advanced Camera for Sur- 
veys (ACS) images in the F435W, F606W, F775W, F850LP 
bands, extensive observations with ESO telescopes both in 
three optical U bands and with J, H, K s near-infrared filters 
and four Spitzer-IRAC images in the 3.6, 4.5, 5.8, 8.0/im 
bands. Morphological parameters (effective radius and Ser- 
sic index) in the F850LP filter and physical parameters (age 
and stellar mass) were already e stimated for each galaxy of 
the sample (|Saracco et al.ll2010l) . Four out of the 11 ETGs 
were already studied in GSL11, thus F606W-F850LP color 
gradients are also available (GSL11). 



3 SURFACE BRIGHTNESS PARAMETERS 
ESTIMATION 

Following GSL11, we have estimated the internal color gra- 
dient of a galaxy as the logarithmic slope of its color pro- 
file, given by the difference between the ^fs50lp{t) and 
A 1 Fi60w(r) surface brightness profiles. As for the F850LP 
and F606W bands, we have modelled the F160W surface 
brightness profile of high-z ETGs with a Sersic law: 

^R)=^ + ^[(r/r.) 1/n -l]. (1) 

To estimate the free parameters of the profile, i.e. the ef- 
fective radius r e (in units of arcsec), the Sersic index n, 
and the normaliz ation term fi e , we have used GALFIT (v2, 
iPeng et al.ll2002r ). This software models the galaxy light dis- 
tribution with a two-dimensional fit. In the fitting proce- 
dure, the galaxy model is convolved with the PSF, provided 
by the user, and the software retains as final solution the 
parameters of the PSF-convolved model that minimizes the 
residuals to the input galaxy image. For each galaxy, we 
have constructed different PSF models, each model being 
obtained by averaging the light profiles of unsaturated stars 
as near as possible to the position of the given galaxy on the 
frame. We have run GALFIT for all different PSFs and have 
retained the case that returns the best residual map (i.e. flat- 
test residuals). For four out of 16 initially selected galaxies, 
we found significant substructures in the residual maps, pos- 
sibly indicating an inaccurate modeling of the galaxy light 
distribution and/or of the PSF. We have excluded these four 
galaxies from the analysis (see Sec. [2}. 

3.1 Robustness of structural parameter estimates 

Given the pixel scale of 0.128" /pixel (0.06" /pixel after driz- 
zling) and the FWHM of the PSF of ~ 0.2", the WFC3 
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Figure 2. Comparison of F160W PSF-convolved Sersic profiles (black curves) for our sample of 11 ETGs, with the observed surface 
brightness profiles measured directly on the given images (red points). Vertical dashed black lines mark the radius of the FWHM of the 
PSF, while red dash-dotted and solid lines correspond to 2 and 3 r e , respectively. In panels related to ETGs 472 and 996, black lines 
correspond to PSF-convolved Sersic profiles derived on the mosaic image of ~ 80ks, while green lines, not visible given the almost perfect 
overlap to the black curves, the same quantity measured on mosaic images of ~ 6ks. 



Table 1. Our sample of galaxies. Column 1: id number Column 2: spectroscopic redshift; Column 3: total magnitude from GALFIT. 
Errors on magnitude are the formal GALFIT errors; Column 4: effective radius in arcsec. A typical error on the estimates of r e is 
0.02 arcsec corresponding to 0.17kpc at z=1.5; Column 5: effective radius in kpc; Column 6: Sersic index; Column 7: F850LP-F160W 
colour gradient. At the median redshift z = 1.5 1 arcsec corresponds to ~ 8.5 kpc. 
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Figure 1. Comparison of effective radius estimated with GAL- 
FfT and 2DPHOT (see the text). The solid line marks the one- 
to-one relation. 



images in the F160W passband are close to the sampling 
limit. Moreover, even with the excellent spatial resolution of 
WFC3, most of the high-z ETGs in our sample have sizes 
comparable to the FWHM. This makes of paramount im- 
portance the accuracy and reliability of the derived PSF 
used to convolve the galaxy model. To assess the indepen- 
dence of our results from the method we have adopted to 
derive the PSF and the consequent robustness of the de- 
rived struc tural parameters, we hav e derived them also with 
2DPHOT (|La Barbera et al.ll200sh . a software working by 
means of different approach than GALFIT. 2DPHOT is a 
fully automatic tool, allowing galaxy surface photometry to 
be performed by fitting galaxy images with two-dimensional, 
PSF convolved, (Sersic) models. Differently from GALFIT, 
2DPHOT constructs a PSF model in a user-independent 
manner, fitting star images on a given frame with a com- 
bination of two-dimensional Moffat functions, taking into 
account PSF asymmetries and image (under-) sampling 0. 
In order to reproduce the spikes of the HST PSF, we have 
modified the 2DPHOT PSF fitting algorithm by smooth- 
ing the (average) residual map from the fitted stars (with 
a 3x3 pixels median smoothing), and adding the smoothed 
residuals to the Moffat-based PSF. The 2DPHOT initial pa- 
rameters for Sersic fitting are also computed through a user- 
independent approach (in contrast to GALFIT), by compar- 
ing the galaxy image with a discrete grid of PSF-convolved 
Sersic models, reducing the problem of spurious conver- 
gences in the (non-linear) optimization procedure. In gen- 
eral, we find a very good agreement between 2DPHOT and 
GALFIT structural parameters. As an example, in Fig[T] we 
compare the effective radii derived with the two softwares. 
Errors in r e has been set equal to 0.02 arcsec, the typical un- 
certainty of this meas ure as we have derived through monte- 
carlo simulations (e.g. lLonghetti et al.ll20 07). No systematic 
difference is found between v e ,GALFiT and r ei 2D phot , de- 
spite of the different approach and methodology adopted. 
Only for one galaxy, with x^^galfit ~ 0.4", the r e estimate 



1 To this effect, at each step of the PSF fitting, the Moffat func- 
tions are convolved with a box kernel, i.e. the pixel of the given 
image. 



differs significantly, with t S: 2dphot being ~ 3 times larger 
than r e ,GALFiT- Given this discrepancy, we have excluded 
this object from our sample thus leaving with 11 galaxies 
in our analysis (see Sec. [2]). To be consistent with GSL11, 
where we estimated structural parameters in the F850LP 
and F606W passbands with GALFIT, we have decided to 
adopt GALFIT estimates of surface brightness parameters 
also for the WFC3/F160W filter. 

To further verify the reliability of surface brightness pa- 
rameters, in Fig. [2] we compare the best-fitting PSF- 
convolved Sersic profiles (black curves), with the observed 
surface brightness profiles (red points), measured in con- 
centric circular coronas of fixed width on the F160W im- 
ages. Errors on surface brightness w ere derived with SEX- 
TRACTOR (|Bertin fc Arnoutsill996D and take into account 
the correlated noise in troduced by the drizzling technique 
IjCasertano et al.l 120001 ). The vertical dashed black lines 
mark the radius of the FWHM of the PSF (0.1 arcsec), while 
red dash-dotted and continuous lines correspond to 2 and 
3 r e , respectively. The figure shows an excellent agreement 
between observed and fitted profiles, from the very central 
region out to a radius of ~ 3r,;, and in many cases well be- 
yond, both for galaxies in the deepest images (ID 472, 996) 
and those in the shallowest fields. This confirms the accuracy 
of PSF modeling and Sersic fitting. For the two galaxies on 
the HUDF area (ID 472 and 996) , the green lines show the 
PSF-convolved Sersic profiles estimated on a mosaic image 
created to match the exposure time of the shallowest images 
available (~ 6 ks; see Sec. [2j. The excellent agreement be- 
tween the black (~ 80 ks) and green (~ 6 ks) solid curves (in 
fact overlapping) in the Figure indicates the homogeneity of 
our measures and proves the feasibility of deriving reliable 
structural parameters even with 2-orbit-long images. 
In Table 1 we report our effective radii, r e (in units of arc- 
sec), and R e (in kpc), as well as the Sersic index and total 
magnitude for each of the 11 ETGs in our sample. 



4 F850LP-F160W COLOR GRADIENTS OF 
HIGH-Z ETGS 

We recall that color gradients are generally defined as the 
logarithmic slope of the color profile: 



FS5QLP-F160W 



A(flF850LP — HF160w)(R) 

A log ft ' 



(2) 



To obtain the measure of this quantity, we have re-estimated 
the surface brightness parameters in the F850LP band, keep- 
ing fixed in the fit of each galaxy the ellipticity and position 
angle to the values derived for the F160W filter. This ensures 
that both the fJ.FS50Lp(r) and /XFi60w( r ) profiles trace the 
variation of light density in the same radial direction. Fol- 
lowing the prescription commonly adopted for nearby ETGs 
and for consistency with GSL11, we have fitted the slope of 
the color profile between 0.1r e ,F850iP and lr e ,F850LP- The 
best-fitting line is obtained by an orthogonal least-squares 
fit, less sensitive to outliers than the direct fit. Fig.|3]reports 
the results. Black curves indicates the color profiles, while 
red lines are the fitted slopes. Dashed lines mark the radius 
of the FWHM of the PSF, while colored area set la errors. 

The values of the gradients and their relative errors are 
reported in Table 1. The error on color gradients is an upper 
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Figure 3. The colour gradients for the galaxies of our sample. Black lines represent the deconvolved colour profiles between 0.1R e and 
R e , and the red lines are the best fitted lines to the models. The dashed vertical lines correspond to the radius of the FWHM of the 
PSF. Colored area indicate the la errors on colour profile slopes. 
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Figure 4. The colour gradients F606W-F850LP and F850LP-F160W (blue and red points, respectively) versus, from left to right, the 
dust extinction Ay, the stellar mass, the redshift, the formation redshift and the global age. 



limit to the true error. We estimated it by associating, at 
each point of the color profile, the error on color, Scol, mea- 
sured directly from the real images in a thin circular corona 
of radius r e (for galaxies with r e >FWHM/2), or FWHM/2 
(for galaxies with r e <FWHM/2). The Scol was estimated 
after reporting the F160W and F850LP images to the same 
pixel scale and PSF (see Sec. 5.1). Our results show that the 
high-z ETGs in our sample have negative or null F850LP- 
F160W color gradients. Indeed, 3 galaxies have color gra- 
dients comparable at la with zero. On the contrary, the 
remaining galaxies show significant negative values which 



reveal the presence of stellar populations redder in the cen- 
ter than those at lR e . The differences observed in the radial 
variation of the F850LP-F160W color within the galaxies 
of our sample show that after 3-4 Gyr at maximum from 
their birth, the stellar content distribution of high-z ETGs 
was not homogeneous. This non homogeneous distribution 
of the stellar content within the galaxies of our sample cor- 
roborates the idea that ETGs have assembled their mass 
through different mechanisms. 

To gain insight into the origin of the observed distribu- 
tion, we firstly have looked for the presence of any cor- 
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relation of observed color gradients with global properties 
of ETGs. Fig. H reports the F606W-F850LP (blue points) 
and F850LP-F160W (red points) color gradients versus the 
dust extinction Ay (first panel), the total stellar mass (sec- 
ond panel), the redshift (third panel) , the f ormation red- 
shift (fourth panel; see ISaracco et al.l (|201lh for the defi- 
nition and estimate of formation redshift) and global age 
(last panel) of the galaxy. No correlation was detected. In- 
deed, the absence of a correlation could be due to the fact 
that we plot color variation within lR e versus quantities 
estimated on the entire galaxy. We have verified that the re- 
sult does not change when the color gradients are estimated 
throughout the whole galaxy. Thus, the various mechanisms 
responsible of the different stellar content observed in our 
high-z ETGs, at the extend of our sample, produce color 
gradients which are independent from redshift, formation 
redshift, stellar mass, global age and dust extinction. This 
result confirms that the observed distribution of color gra- 
dients is primarily due to differences in the stell ar content 
of hig h-z ETGs. Different results were reached bv lGuo et al.l 
l)201ll ) who found a correlation between dust and color gra- 
dients in the direction o f steeper gradien t for galaxies with 
higher dust-obscuration. iGuo et al.l (|201ll ) fixed the dust ex- 
tinction assuming a Salpeter IMF in the SED fitting, while 
we have adopted a Chabrier IMF. Furthermore they derived 
their color gradients fitting the color profiles between ~ lr e 
up to 5-8r e . To properly compare our results with those of 
Guo et al., we have re-estimated the color gradients in our 
sample in regions similar to those selected in their paper, 
and computed Ay assuming a Salpeter IMF (see Table 2). 
Even in this case we do not find any correlation between 
dust extinction and color gradients. Actually, two galaxies 
of our sample were analyzed also in Guo's paper (ID 472 and 
996 our paper, 23555 and 22704 in their). Our estimates of 
the effective radius of these two galaxies are ~ 13% smaller 
than their ones, thus consistent within the typical errors on 
this quantity. Moreover, for both galaxie s our estimates o f 
Ay a re consistent with those reported bv lGuoet aDl|201ll ). 
Since IGuo et all l|201ll ) did not report the gradients values, 
we try to qualitatively compare our color profiles with their 
ones. We would like to stress here that, for these galaxies, 
the typical error on color at 7r e is ~ 0.15 mag. For the 
galaxy 472 (see our Fig[7] and magenta line in the lower 
panel of their fig. 5) the agreement between the two color 
profiles seems very good. For the galaxy 996 we have found 
a global negative trend in color values, but the amplitude 
of the color gradient we detect is less steep than the one 
measured by Guo, due to the external region (5r e <r< 7r e ) 
where their profile falls suddenly down, while our one follows 
a quite costant variation (see our Fig[7]and green line in the 
lower panel of their fig. 5). The reason for this difference 
can be due to the presence of another galaxy at very low 
projected distance (see our Section 5.1 and their Fig. 1). In 
the F850LP band, this source is not detected due to its faint 
magnitude, thus the extreme blue color observed in the ex- 
ternal region by Guo et al. could be due to an underestimate 
of the flux in the solely F160W band. On the contrary, in the 
fitting procedure of the surface brightness profile, we have 
simultaneously modeled both galaxies and assured from the 
residual map that the contribution of the other galaxy was 
properly taken into account. If this explanation is correct, 
the consequent flattening of the color gradient value of this 



galaxy would strongly weaken the correlation between color 
gradients and extinction also in their sample and the even- 
tual residual correlation should be totally due to the small 
statistics of their sample. 



5 SPATIAL DISTRIBUTION OF STELLAR 
POPULATIONS IN HIGH-Z ETGS AND 
RADIAL VARIATION OF THEIR 
PROPERTIES. 

In order to constrain the possible physical mechanisms re- 
sponsible of the mass accretion within the ETGs, we have 
designed a new method able to investigate the spatial distri- 
bution of the underlying stellar populations in high-z ETGs 
which is directly correlated to the mass assembly processes. 
Briefly, we have derived the color maps of our galaxies and, 
on their basis, we have defined for each ETG of our sample 
an internal and an external region. Then, we have modelled 
the global stellar content of each ETG as formed by two 
main stellar components, one dominating the internal re- 
gions, and the other the external ones. With synthetic mod- 
els of composite stellar populations we have constrained the 
stellar properties of the two components in order to simul- 
taneously reproduce the observed color gradient(s) (in one 
or two colors) and the whole global SED of the galaxy from 
0.3/im to 8/im. 

5.1 Color maps of high-z ETGs: a direct look at 
the spatial distribution of their stellar 
content. 

To have a direct look on the punctual 2-dimensional dis- 
tribution of F850LP-F160W color in high-z ETGs, we have 
derived their color maps. To this aim, we have ri-reduced 
the WFC3 images to a pixel scale of 0.03 "/px to match 
the one of F850LP-mosaic images, and used the IRAF tasks 
geomap/geotran to align the two mosaics. Then, we have 
degraded the F850LP images (FWHM ~ 0.12") to the same 
PSF of the F160W images (FWHM ~ 0.2"). To this aim, 
we needed to derive the kernel function K(r) which regulates 
the transformation between the PSFs in the two bands, i.e. 
which holds: 

PSF FS soLp(r) * K{r) = PSF F160W (r) (3) 

where PSFfssolpM and PSFFi60w(r) are the point spread 
functions in the F850LP and F160W band, respectively, and 
where the symbol ★ denotes a convolution. For each galaxy, 
we have modeled its PSFfssolpW and PSF_Fi6iw(r) func- 
tions and through Eq[3]we have recovered its kernel function. 
The PSFF850ip(r) and PSFFi60w(r) functions of a given 
galaxy have been derived by fitting the light profiles of the 
unsaturated stars we have adopted as PSF in GALFIT and 
that we have already tested to be a good approximation of 
the true PSF for that galaxy. In Figure [5] we report a repre- 
sentative example of the fractional encircled energy for the 
real PSF (black line) of a star in the F160W (left panel) and 
F850LP (central panel) band and the corresponding fitted 
model (red line). The agreement between the two curves is 
extremely good in both cases, confirming the quality of our 
fit. In the last panel of Fig. [5] we report the encircled en- 
ergy distribution of the PSFFi60w(r) (black line) and of the 
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Figure 5. Left panel: the fractional encircled energy for the real PSF in the F160W image (black line) and for the model (red line). 
Central panel: the same for the F850LP image. Right panel: the fractional encircled energy for the real PSF in the F160W image and 
for the PSF of the F850LP image convolved with the kernel function K(r). 




Figure 6. F850LP-F160W colour maps for the galaxies of our sample. The maps are color coded, indicating, from blue to red, region 
of the galaxy characterized by a redder emission. For each galaxy, the inner white circle constrain the internal region, while the two 
external ones define the circular corona which identify the external region. Magenta, blue and red circles set the effective radius in the 
F606W (when available), F850LP and F160W band, respectively. The white line corresponds to 1 arcsec. 



PSF_F850Lp(r) model convolved with the kernel (red line). 
The match between the two functions is perfect and differ- 
ences < 2% appear at R> 4", when present. In particular, 
it is to note the optimal agreement in the central region 
that, as said before, is the one mainly affecting the galaxy 
light profile. This has allowed us to degrade the F850LP im- 
ages to the same PSF of F160W images and to derive the 
color maps presented in Fig. [B] The images are color coded 
with the reddest color indicating the region of the galaxies 



with highest F850LP-F160W values. Red, blue and magenta 
circles indicate the r e in the F160W, F850LP and F606W 
band, respectively. As pointed out before, we have at dis- 
posal information on F606W band only in 4 galaxies out of 
11. In some cases the r e in the F160W band is represented 
with a black circle being the red one not visible on the color 
map. The white lines correspond to 1 arcsec. In Fig. [6] we 
have marked all the pixels 4.5cr above the mean background 
in the F160W images as belonging to the galaxy. In spite of 
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this high threshold, the color maps extend to regions well 
beyond 2r Gi F850iP and at radius greater than 4r 6i F850LP for 
5 cases (12294, 996, 2361, 2111, 472). The color map of the 
whole galaxy show that differently from the gradients es- 
timated up to lr ej F850i,p, oil high-z ETGs of our sample 
have a negative color gradient. Moreover, they point out 
that although the general trend is the same, the F850LP- 
F160W color varies from the internal to the external regions 
with different gradients showing that both the spatial dis- 
tribution and the properties of the stellar content are not 
homogeneous in high-z ETGs. 

5.2 Radial variation of the stellar populations 
properties of high-z ETGs. 

To identify the stellar population parameters whose varia- 
tion can account for the observed color gradients, we have 
used the following procedure. On the basis of the color maps 
which show a quite radial variation of color within the ETGs 
of our sample, we have schematized the whole stellar con- 
tent of each of them with a two-components model assum- 
ing that a component dominates the central regions and the 
other dominates the external regions. By means of the color 
maps we have selected the two areas (white circles in Fig. [6]) 
in order to maximize both the difference in the mean color 
of the two regions and the radial distance among them. At 
fixed gradient, maximizing the distance of the two regions 
means maximizing their color difference, and this will facil- 
itate in detecting differences between the stellar population 
parameters of the two components and hence in detecting 
possible radial variations. 

We have assumed as initial guess for the stellar param- 
eters of the two components the values derived from the 
fit of the global SED [Age 9 ; ofjai , Z globa i, A v , globa i, T g i oba i]. 
Briefly, for each galaxy, the global stellar population param- 
eters were derived through the fit of the spectral energy dis- 
tribution defined by 14 photometric points and with known 
spectroscopic reds hift. The fit was perfor med with the soft- 
ware Hyperzmass ( Bolzonclla et all 2000) using the compos- 
ite stellar population models of Chariot & Bruzual (here- 
after CB07) with an exponentially declining star-formation 
history (oc exp ~ t ^ T ), and a Chabrier initial mass function 
(|Chabrierl 2003) The best solution was selected on a grid of 
model with varying age, star-formation time scale r [0.1, 0.3, 
0.4, 0.6] Gyr, and solar metallicity. Extinction Ay was left 
as free parameter in the 0.0 - 0.6 range and extinction curve 
of ICalzetti et all (12 000) was adopt ed. For more details on 
SED fitting see lsaracco et alj ()2010h . For each galaxy of our 
sample, the set of the four best-fit parameters are reported 
in Table [2] 

In order to reproduce the color gradient we observe, it 
is necessary that one or more parameters of the two compo- 
nents vary from the internal to the external region. Unfortu- 
nately, the effect of age, metallicity, dust and r variation on 
the galaxy emission is degenerate in the spectral region we 
are observing, preventing, in fact, to identify the contribu- 
tion of each parameter to the total emission. For this reason, 
we have performed our analysis by investigating the effect of 
each single parameter in turn, without considering the pos- 
sibility of either a correlated or an anticorrelated variation 
of age and metallicity within galaxies. Thus, we have fixed 
in both components three out of the four parameters in the 



fitting. For example to investigate radial age variation as 
possible driver of color variations, we fixed the metallicity, 
the dust extinction and the star-formation time scale to the 
value obtained by the fit of the global SED: Z g i oba i, Ay, global, 
Tgiobai- We have looked for the value of age to be associated 
to the internal and to the external populations, Age™ and 
Age ou t respectively, which best reproduce the observed color 
gradients. In particular we have chosen the value of Age™ 
so that minimizes the quantity 

(F850LP - Fl&QW) mod ,in - (F850LP - F160W) obs>in \+ 
+ (F606W - F850LP) roO(M n - (F606W - F850LP) obs ,in\ 

(4) 

where (F850LP-F160W) mod , m and (F606W-F850LP) mod , in 
are the F850LP-F160W and the F606W-F850LP colors pre- 
dicted by the model defined by the stellar population pa- 
rameters (Age in , Z g i bai, ^v, a iobai, T globa i) , and (F850LP- 
F160W) o i, s , in and (F606W-F850LP) oi)s , in are the mean 
value of the colors observed in the internal region. Similar 
procedure has been adopted to identify Age out . In the cases 
for which F606W-F850LP color gradients are not available 
we clearly ignore the second lines in Eq(4] To fix the con- 
tribution to the total stellar mass of the two components so 
selected, we fit the global SED sampled by 14 photometric 
points with a linear combination of the SEDs of the two 
populations. 

We have repeated this analysis also for variation of 
metallicity Z, keeping fixed [Age g i b a i, Av, g ioba,l, T g i oba i], and 
star-formation r with [Age g i oba i, Z g i oba i, A^j| oW ] fixed. We 
have looked for the values of metallicity to be associated to 
the internal and external population, Z; n and Z out respec- 
tively, on a grid of sub-solar and super-solar values: 0.2Zq, 
O.4Z0, Zq and 2Zq, while for the star-formation time scale 
we have considered the range 0.1 -0.6Gyr. 



5.2.1 An example 

In Fig. [7] we report an example of the analysis described 
above for a galaxy of our sample. We select a case which 
presents some peculiarities to guide the interpretation of the 
analysis in this more complex case and left to the reader the 
interpretation of the analysis of the other galaxies not dis- 
cussed here but presented in electronic form. The first col- 
umn reports the analysis relative to the effect on galaxy 
colors of pure age variation from the inner to the outer 
regions. The second and the third column show the same 
but for metallicity and r radial variation, respectively. In 
each column we report, from top to bottom, the F850LP- 
F160W color gradient, the F606W-F850LP color gradient, 
when available, and the global SED. In the panels related to 
the color gradients, the black lines are the color profiles. Red 
lines are the fit to the color profiles up to the external regions 
derived following the same method describe in Section 4. 
Red points are the mean color values of the internal regions 
measured through color profiles, i.e. (F850LP-F160W) o t s ,i„ 
and (F606W-F850LP) oi)jMn in Eq. 0] Blue points are the 
same for the external regions. The horizontal error bars in- 
dicate the extension of the two areas we select (white circles 
in color maps) and over which we computed the mean colors. 

In the bottom panel we report the observed SED (black 
dots), the best fitting template (black line) and the relative 
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Table 2. Stellar population parameters derived from the fit of global the SED assuming a Chabrier IMF and a Salpeter IMF. Column 
1: id number; Column 2: Age obtained with a Chabrier IMF|Age obtained with a Salpeter IMF; Column 3: r obtained with a Chabrier 
IMF|r obtained with a Salpeter IMF; Column 4: Ay obtained with a Chabrier IMF|Av obtained with a Salpeter IMF; Column 5: log 
M + obtained with a Chabrier IMF| logM* obtained with a Salpeter IMF. In each fit we assume solar metallicity. 
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best fitting parameters (black text). Red and blue lines are 
the contribution of the internal and external population, re- 
spectively, to the total emission, and hence to the total stel- 
lar mass. The amount of their contribution was fixed such 
that their linear combination (magenta, cyan or green lines 
for age, metallicity and r variations, respectively) minimizes 
the residual with the observed SED. In each panel we report 
the x 2 value of this fit. In red we present the value of the 
varying parameter which minimize the Eq[4] as well as their 
contribute in mass, and in blue the same for the external 
region. 

The solely linear radial variation of this parameter from 
the internal value (red text), to the external value (blue 
text) will produce a F850LP-F160W and a F606W-F850LP 
color gradient indicated in the top and middle panels by the 
dashed lines (magenta, cyan and green in the case of age, 
metallicity and r variation) Fig. [7] shows that the color gra- 
dient we observe in this galaxy can be due to a variation of 
the age of the stellar populations from the inner to the outer 
regions. Indeed, while the global fit returns us a mean age 
of 1.9 Gyr (black text in the bottom left panel), our analysis 
shows that a population of 2.7 Gyr dominating the central 
region (red text in the bottom left panel) and a younger pop- 
ulation of 1.0 Gyr dominating the external regions (blue text 
in the bottom left panel), can perfectly reproduce within the 
errors (magenta dashed line in top and middle panels) the in- 
ternal and external observed F606W-F850LP and F850LP- 
F160W colors, and hence the color gradients. Further, the 
analysis of the SED suggests that in this galaxy the older 
stellar population contributes to the stellar mass for more 
than 90% (red line, bottom left panel). On the contrary, 
neither a pure metallicity variation, nor a variation of the 
star-formation time scale alone can reproduce the observed 
color gradients. In the case of metallicity (middle column), 
the internal colors are well reproduced by a population with 
solar metallicity (red text in the bottom central panel), i.e. 
the same of the global fit, while the external color cannot 
be simultaneously reproduced by a single value of Z. 
For what concerns the analysis of the r variation this galaxy 
requires particular attention. The fit of the global SED sug- 



gests a mean r for this galaxy of O.lGyr, that is the lower 
value in the grid of the star-formation time scales consid- 
ered. To cope with this, we assume a population of 2.7 Gyr 
with t = 0.1 Gyr as representative of the inner part, as we 
already found that this parameter well fits the inner colors. 
For the external regions, fixing the age to 2.4 Gyr, the value 
of t that best reproduces both external colors is 0.4 Gyr. 
However, despite it is the best fitting value, it is not able to 
reproduce within la the observed F850LP-F160W external 
color. 
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Figure 7. First column: analysis of radial age variation as possible driver of observed color gradients. Top and middle panels show the 
F850LP-F160W and F606W-F850LP color profiles, respectively (black lines), whit the relative fit (red lines). The bottom panel shows 
the observed global SED (black dots) as well as the stellar population parameters (black text) derived from its fit (black line). Among all 
the stellar populations with varying age but star-formation time scale, metallicity and dust extinction equal to the value obtained from 
the fit of the global SED, the red text identifies the one that best simultaneously reproduces the colors observed in the internal region 
(red dots in the top and middle panel). Similarly, blue text reports the value of the age of the stellar population that best reproduce the 
external colors (blue dots in the upper panels). The pure linear variation of the age of the galaxy stellar populations from the internal 
value to the external one, will produce the F850LP-F160W and F606W-F850LP color gradients represented with the dashed magenta 
lines. The contribution to the total stellar mass of the two populations is fixed in order to best fit the whole SED and are reported in 
the bottom panel (red and blue lines). Second and third column: the same of first column but for metallicity and star-formation time 
scale variation. 



Table 3. Color gradients produced by the variation of a single stellar population parameters from the internal to the external regions. Column 1: top line: galaxy id number; bottom to 
line: compactness C defined as R e ,z=o /Re where R e is the effective radius of the galaxy and R e , z= o is the radius that a galaxy of equal stellar mass would have at z=0 as derived by the 
size mass relation of lShen et al.l d2003T) - ETGs with effective radius more than one sigma smaller than those predicted by the local size-mass relation, i.e. our compact galaxies, turn out ^ 
to have C^2. ; Column 2: Radius up to which we extend the analysis (R ex t, hi R e F850LP unit); Column 3: F850LP-F160W measured color gradient |F606W-F850LP measured color 
gradient derived up to K ex t from the fit of the color profile; Column 4- top line: F850LP-F160W color gradient |F606W-F850LP color gradient produced by a solely radial age variation 
from Age; n in the internal region to Age otl < in the external region (see bottom line); bottom line: Agei n : age of the internal population - Age ou t: age of the external population - V age : © 
age gradient derived as cilog Age/dlog r; Column 5: the same of column 4 for metallicity radial variation; metallicity gradient is derived as cflog Z/dlog r; Column 6: the same of column 
4 for star-formation time scale radial variation; Column 7: the same of column 4 for IMF's slope radial variation. In the case of star-formation time scale and IMF's slope radial variation S_ 
the V T and V a are omitted. In bold case are reported the variations able to reproduce both the color gradients and observed global SED. ° 
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-0.86| — 
3.0 - 0.8 - -0.84 


-0.82| — 
2Z Q - 0.2Z Q - -1.47 


-0.52| — 
0.1 - 0.6 


0.07| — 
1.5 - 3.5 


2543 
2.6 


2.0 


-0.73|— 


-0.68 — 
3.5 - 1.8 - -0.33 


-0.67| — 
2Z Q - 0.2Z Q - -1.19 


-0.42| — 
0.1 - 0.6 


-0.41 — 
3.5 - 1.5 


2111 
3.8 


4.9 


-0.76|0.05 


0.17|0.3 
1.0 - 1.1 - +0.05 


0|0 

Zq - Z Q - 


0|0 
0.1 - 0.1 


-0.01j0.02 
3.0 - 3.5 


472 
5.4 


8.6 


-0.31|0.83 


0.22|0.93 
0.9 - 1.4 - 0.18 


0.21|0.71 
Zq - 2Z - 0.28 


0|0 
0.1 - 0.1 


-0.08|0.30 
1.5 - 3.5 
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The fact that we have changed the global age value 
does not influence our analysis. Indeed, we make the starting 
assumption to associate to the two components the stellar 
population parameters of the global fit since it can be a 
good representation of the whole galaxy. However what we 
are interested in is the variation of these parameters as 
possible driver of the observed color gradients, and not their 
absolute values. In Table[3] we report the summary of the 
results obtained for the other galaxies of the sample, while 
the plots related to their analysis similar to Fig[7]and with 
the same conventions are available in electronic form. 

5.2.2 Dust distribution and color gradient 

Actually, as pointed out before, one of the causes of the 
radial variation of color we observe in our ETGs can be a non 
homogeneous distribution of the dust within the galaxy. Fig. 
|4]shows that for half galaxy of our sample the fit of the global 
SED returns a value of Av = 0. Thus, for these galaxies (ID 
23, 11888, 2361, 2196, 472) seems not plausible a scenario 
whereby the main driver of the color gradients is a radial 
distribution of dust. For the remaining galaxies, the observed 
global SED was best fitted assuming the presence of a non 
null quantity of dust. As in the case of age, metallicity and r 
we test if, starting from a flat color distribution, a pure radial 
variation of the dust extinction (more dust in the center than 
in the outskirt) can reproduce the observed color gradients. 

In particular, we have made the working hypothesis that 
the ETGs of the sample have no dust in the outskirt. Thus, 
the colors we measure in the external regions (i.e. (F850LP- 
F160W) O 6s, out) are the intrinsic colors of the galaxies. In 
these hypothesis, the reddest colors observed in the center 
of the galaxies (i.e. (F850LP-F160W) o b s , in ) are the direct 
effect of the extinction of the dust. The necessary amount 
of dust extinction to redden the intrinsic colors to those 
observed in the center is: 

A v = 

[(F850LP - FimWUsM - (F850LP - F160W) obs , out ] x R v 
K(F850LP) - K {FimW) 

(5) 

where K(F850LP)A V /R V and K(F160W)A V /Rv are the 
extinction in both bands due to the dust. We have derived 
K(F850LP) and K(F160W) following ICalzetti et all (|2000H 
and have assumed Rv = 4.05 ± 0.80. For the galaxy 2111, 
Eq. \5\ leads to a dust extinction Av in the center ~ 1. This 
value is a lower limit. Indeed, if dust was not localized only 
in the center as we have assumed, but also in the periph- 
ery, the amount of dust necessary to generate the observed 
color gradient should be even higher. The fit of the global 
SED returns a global value of Av = 0.35, much lower than 
the one we obtain. Moreover, the amount of dust predicted, 
that necessarily affects also the emission in the F606W and 
F850LP bands, produce F606W-F850LP color in the central 
region of the galaxy redder than the observed one. Thus, for 
this galaxy, dust does not seem the main driver of the ob- 
served color gradients. 

For galaxies 2286, 2239, and 2543 the analysis indicates 
that, in our hypothesis, their color gradients can be obtained 
assuming a dust extinction in the central regions of Av = 
0.55, 0.29 and 0.83, respectively. Table [2] shows that for the 
galaxy 2239 the dust extinction we need to reproduce the 




Mass [MO] 

Figure 8. The relative amount of low-to-high mass stars pre- 
dicted by 5 different IMFs with analytic form dN/dM a M" a 
and a = 1.5, 2.0, 2.35, 3.0, 3.5. 



color gradient is lower than the one predicted by the fit of 
the global SED. On the contrary, galaxies 2286, 2543 require 
a dust extinction ~ 1.5-2 times higher than the one we ob- 
tained from the global SED fitting. These values confirm 
that for these galaxies, if not the main driver, a significant 
contribution of dust in generating the observed color cannot 
be excluded. 

Finally, for the galaxy 996 and 12294 we have obtained 
Ay = 1.08 and 1.53 respectively, results which disagree with 
the best-fit values of 0.10 and 0.25. In these cases the dust 
extinction we need to reproduce the observed color gradients 
is 6-10 times higher than the one obtained from the global 
SED fitting, suggesting that, although it is not possible to 
exclude it at all, the radial variation of dust extinction is 
certainly not the main driver of the color variation we ob- 
serve. 

5.2.3 Radial variation of initial mass function as possible 
driver of color gradients 

Among the main hypothesis to explain the color gradients 
we observe, there is also a change in the initial mass function 
(IMF) from the inner to the outer regions of the galaxy. 

In spite of the fact that a radial variation of IMF could 
be a reliable way to produce color gradient, it has never been 
explored due to the lack of suitable stellar population models 
For this reason a grid of stellar population models has been 
kindly produced ad hoc by Stephan Chariot for this purpose: 
reproducing the spectral energy distribution of galaxies with 
different proportion of low-to-high mass stars. In particular 
we have assumed an IMF analytically described by a power 
law form: 

dN/dM oc M~ a (6) 

with 5 different values of a: 1.5, 2.0, 2.35, 3. 0, 3.5. The valu e 
2.35 corresponds to a pure Salpeter IMF l|Salpeterl Il955l ). 
Fig. [8] shows the relative different amount of low-to-high 
mass stars produced by the 5 IMFs. 
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The evolution of a star with time is strictly related to 
its mass. Thus, at fixed age, metallicity, star-formation time 
scales, and dust extinction, the spectral energy emission of a 
galaxy will change with the IMF due to the different amount 
of low-to- high mass stars. Left panel of Fig.[9]shows how the 
template of the global SED of a possible galaxy of our sam- 
ple (age = 3.5 Gyr, r = 0.5 Gyr, z=1.6) would change with 
different IMFs. To observe the relative differences, we have 
normalized the observed SEDs in the wavelength range [8250 
A- 8750 A]. The colors of the different SED follow the same 
code of Fig. [8] Since we normalized all the curves approxi- 
mately in correspondence of the F850LP filter, when looking 
at the emission around 16000 A this plot returns immedi- 
ately how the F850LP-FI60W color changes with the IMF. 
In particular, for this template, IMF with an higher abun- 
dance of low mass stars (red line) returns redder color than 
the one with lower number of low-mass stars. Actually, the 
color of a star changes through all its life cycle, as it moves 
on the HR diagram according to its mass. So, the contri- 
bution of low and high mass stars to the total emission in 
a fixed band changes at the variation of the ratio age/r. 
Right panel of Fig. [5] shows the same plot on the left but for 
a galaxy with a different ratio Age/r. This simulated galaxy 
has age = 1.1 Gyr and r = 0.1 Gyr. Even in this case, we 
normalize the curves in the wavelength region [8250 A- 8750 
A]. This plot shows how in this case, the near-IR domain 
is dominated by the light coming from the high-mass stars 
contrary to the previous case. These plots clearly emphasize 
two main aspects: colors are affected by the IMF, hence the 
radial color variation we detect can be effectively due to ra- 
dial variation of IMF, and the shape of the SED depends on 
the ratio age/r and IMF. 

To investigate the hypothesis of an IMF radial variation as 
the main driver of color gradient, we have adopted the same 
method used to study age, metallicity and r radial varia- 
tions. For each galaxy we have fixed the stellar population 
parameters of the internal and external populations to those 
derived from the global fit and have looked for the value of 
the IMF's slope able to reproduce the observed color gra- 



dients. Since our goal is to investigate the radial variation 
of low-to-high mass stars abundance, we have re-fitted the 
observed global SED assuming a Salpeter IMF instead of 
the Chabrier. In Tabled the result of the SED fitting with 
a Salpeter IMF are reported. We have assumed the inter- 
nal and external populations defined by these new global 

parameters (A g i bal,Sal, global, Sal, Av,global,Sal, T g tobat,Sah 

a=2.35) and we have looked for the value of a to be associ- 
ated to the internal and external population, atin and a out 
so that the two populations defined by the set of param- 
eters [Aglobal,Sal, ^global, Sal, Ay,glabal,Sal, T global, Sal, Clin) 
and (A g lobal,Sal, 7iglobal,Sal, Av,global ,Sal , T global, Sal, Clout) 

best reproduce the observed color gradients and the global 
SED. 

In all but two cases a variation of the abundance of low- 
to-high mass stars can be excluded as the main driver of 
the observed color gradients. On the contrary Fig. HOI shows 
that for the galaxy 2239 an external population account- 
ing for the great part of the mass (> 90%) formed with an 
IMF with a = 3.0, i.e. with an higher number of low mass 
stars than a Salpeter, and a small contribution in mass of 
an internal population with an IMF defined by an a = 1.5 
can simultaneously reproduce the F850LP-F160W observed 
color gradient we have derived and the observed global SED. 
At the same way, the F850LP-F160W color gradient and 
the global SED of the galaxy 2543 are well reproduced by 
a very massive population dominating the internal region 
characterized by an IMF with a — 3.5 and by a small con- 
tribution of an external population described by an IMF 
with a — 1.5. As noted before, the different ratio age/r of 
these two galaxies reflects in inverted radial trend of IMF's 
slope to reproduce negative color gradients. Thus, for these 
two galaxies a variation of the slope of the IMF, and hence a 
variation in the abundance of low-to-high mass stars can be 
a possible driver of the radial color variations we detect. In 
Tabled the results of the analysis for all the galaxies of the 
sample are shown, while in electronic form we present their 
plots. 
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the global SED. For four of these galaxies (ID 2286, 2239, 
2196, 2543), also radial metallicity variation can reproduce 
the color gradient. On the contrary, a radial variation of the 
star-formation time scale, dust content, and abundance of 
low-to-high mass stars are able to reproduce both color gra- 
dients and global SED in only few of these 6 galaxies. More- 
over, we pointed out that a radial variation of the slope of 
the IMF is not able to reproduce the observed radial color 
variation and spectral emission in any of the four galaxies 
for which both F606W-F850LP and F850LP-F160W color 
gradients are available. 

These results assign to age and metallicity a central role 
in generating the observed color gradients. As pointed out 
before, a contribution of other parameters to the color vari- 
ations we detected cannot be ruled out. 
For the remaining five galaxies (ID 2111, 11888, 12294, 472, 
2361) the variation of a single stellar population parame- 
ter is not able to reproduce the observed colo r gradients 
and g lobal SED consistently with the findings of lGuo et al.l 
l|2011 ) on their sample. This suggests that for these galaxies 
a simultaneous variation of several parameters has to be in- 
voked to reproduce the observed color variations and global 
SED. For these galaxies, an approach that investigates the 
simultaneous variation of more than one parameter should 
be considered. 



Figure 10. Analysis of the radial variation of IMF's slope as 
possible driver of observed color gradients. Top panel shows the 
F850LP-F160W color profile (black line) whit the relative fit (red 
line). The bottom panel shows the observed global SED (black 
dots) as well as the stellar population parameters (black text) 
derived from its fit (black line) assuming a Salpeter IMF. Among 
all the stellar populations with varying IMF's slope but age, star- 
formation time scale, metallicity and dust extinction equal to the 
value obtained from the fit of the global SED, the red text identi- 
fies the one that best reproduces the color observed in the internal 
region (red dot in the top panel). Similarly, blue text reports the 
value of the slope of the IMF that best reproduce the external 
color (blue dot in the upper panel). The pure linear variation of 
the IMF's slope from the internal value to the external one, will 
produce the F850LP-F160W color gradient represented with the 
dashed yellow line. The contribution to the total stellar mass of 
the two populations is fixed in order to best fit the whole SED 
and are reported in the bottom panel (red and blue lines). 



6 SUMMARY OF THE RESULTS 

We have found that ~ 70% (8 out of 11) of the ETGs in 
our sample show negative color gradients at more than ~ 2a 
(>4<t in 5 ETGs) in the 0.1R e -lR e range, the effective ra- 
dius range usually adopted to study color gradients in local 
ETGs. The remaining 30% show a gradient consistent with 
a flat color profile. Extending the analysis at R>R e , enclos- 
ing the whole galaxy, we have found that the fraction of 
high-z ETGs with negative F850LP-F160W color gradients 
rises up to 100%. In fact, we have generally found a steep- 
ening of color gradients extending the fit to regions out of 

Re (2R e <R< 8-Re) 

For 6 galaxies of our sample (ID 23, 2286, 2239, 2543, 2196, 
996) a solely radial variation of the age of the stellar pop- 
ulations can simultaneously reproduce at less than la the 
observed color gradients throughout the whole galaxy and 



7 DISCUSSION AND CONCLUSIONS 

The results of our analysis establish the heterogeneity of the 
stellar content in high-z ETGs. For ~ 50% of the galaxies 
of our sample (6 ETGs) we have found that a pure radial 
age variation, with the oldest stars located at the center of 
the galaxy, is able to reproduce the observed radial color 
profile and global SED. The age gradients we have detected 
for these galaxies span a range from -0.84 to -0.28 dex per 
radial decade. For 4 out of these 6 ETGs the color gradients 
we have observed can be also accounted for by solely metal- 
licity gradients whose strength varies from -0.35 to -1.45 dex 
per radial decade. 

Due to the lack of similar measurements on other sample of 
high-z ETGs we have compared the age/metallicity gradient 
values we have found with those observed in the local ETGs. 
Studies of ETGs at lower redshift affirm that color gradients 
are mainly due to radial metallicity variations. Assuming 
the age of the stellar populations constant throughout the 
galaxy, as we did in the case of metallicity gradient anal- 
ysis, color gradients in local ETGs turn out to be repro- 
duce d by a mean metallicity gradients ranging from -0.16 -f- 
-0.3 llPeletier et al.lll990l ; lldiart et al]|2003l ; lTamura fc Ohtal 
2003). A further confirmation of the metallicity as main 
driver of local gradients comes out from studies at interme- 
diate redshift which show that the color gradients evolution 
is better accoun ted for by the pa s sive evolution of m etal- 
licity gradients (|Saglia et al.ll2000l ; Irlinklev fc IrrJl200ll ). In 
fact, our results seem to be inconsistent with these find- 
ings. In our sample of high-z ETGs we have found that only 
4 galaxies out 11 have color gradients well reproduced by 
pure radial metallicity variation. Moreover, the metallicity 
gradients we have detected in high-z ETGs ([-0. 3-1— 1.45]) are 
systematically steeper than tho se typically observed in lo- 
cal ETGs ([-0.16-7-0.3]), even if lOgando etafl (|2005h found 



16 A. Gargiulo et al. 



that this range become wider ([0.0 -1.0]) fo r nearby mas- 
sive (M* > 1O 1O M ) ETGs see also (see also ISpolaor et all 
2009). The steeper metallicity gradients that we have de- 
tected derive by a radial metallicity variation from superso- 
lar (2Zq) values in the inner regions to subsolar values in the 
external regions (0.2Zq). Although such ext reme values of 
Z have been observed als o in few local ETGs l|Mehlert et al.l 
120031 ; lRicke7 et al.l [2008) . our results show that metallicity 
gradients in high-z ETGs of our sample are only marginally 
comparable with the typical metallicity gradients detected 
in local ETGs. This result seems to point in favor of a pos- 
sible evolution of the metallicity gradient in the last 9Gyr. 

Concurrently, studies on cluster ETGs at low and in- 
termediate redshift show that pure age variations in their 
stellar po pulations are not able to account for their color 
gradients (|Saglia et al.|l2000l ). In contrast to local results, we 
have found that for ~ 50% of the galaxies of our sample a ra- 
dial variation of stellar populations age alone can reproduce 
the observed color gredients and global SED. In fact, recent 
studies investigating the simultaneous radial variation of 
both age and metallicity confirm metallicity gradients (~-0.4 
dex per radial decade) as the main driver of observed color 
gradients in local ETGs but found also a small contribution 
to color variati on of positive age gradient (^0 . 1 dex per 
radia l decade) (|La Barbera fc de Carvalhol 120091 ; IWu et all 
2005). The age/metallicty degeneracy affecting optical col- 
ors does not allow us to consider the simultaneous radial 
variation of both parameters and hence to detect a possible 
positive age gradient, whose presence in local ETGs, actu- 
ally, is still matter of debate. 

To spread light on this issue, high-z ETGs constitute the 
ideal place to investigate the presence of an age gradient. 
Indeed, at fixed radial variation of age Aage, its effect on 
color profile is much more enhanced when stellar populations 
are younger, hence in high-z ETGs. This effect is clearly 
shown in left panel of Fig. 1111 where we report the differ- 
ences observed in the F850LP-F160W color of two stellar 
populations with age which differ of 2Gyr (r = O.lGyr, black 
solid curve), as a function of the age of the youngest stel- 
lar population. The same Aage produces a difference in the 
F850LP-F160W color of the two populations that is ~ 10 
times higher if observed in high-z ETGs (Age < ~ 4Gyr) 
respect those observed in local ETGs (Age ~ lOGyr). A 
typical radial variation of 2 Gyr, as the one we measure in 
the ETGs of our sample, will produce in a local ETG a vari- 
ation in the F850LP-F160W color of ~ 0.05 mag, thus at 
the very limit of the detection. On the contrary, the same 
age variation will result in a color variation of ~ [0.3-0.5] 
mag for high-z ETGs. Red and blue lines report the same of 
black line, but for pure metallicity variations. In particular, 
red line show the variation in the F850LP-F160W color ob- 
served in two populations with metallicities 2Zq and O.2Z0, 
while blue line in two populations with metallicities Zq and 
O.2Z0. In the right panel of Fig. 1111 the color gradient that 
the above age/metallicity variations would produce when 
occurring between 0.1R e and 3R e . Differently from age vari- 
ation, the effect of a metallicity variation on color, and hence 
on gradient, increases with the age of the galaxy of a factor 
~ 2 from high-z ETGs to local ETGs. These plots emphasize 
how challenging is the detection of an age gradient in local 
ETGs due to its almost negligible effect on color profile. On 
the contrary, in high-z ETGs age and metallicity variations 



produce comparable effect on color profile, thus setting the 
ideal condition for their detection. This comparison with lo- 
cal samples is only meant to have an indication on how the 
results obtained at high-z, both in terms of age/Z gradients 
and in terms of internal and external age/Z values, relate 
with the local values. On the other hand, the unknown evo- 
lution experienced by ETGs in ~9 Gyr from z~1.5 to z=0 
can affect the stellar properties and distribution (e.g. minor 
mergers triggering secondary burst of star formation) mak- 
ing complex any comparison with local universe. In fact, to 
properly face on high-z and local values samples of ETGs 
selected in an homogeneous way, a dataset able to trace the 
evolution of the same rest-frame color gradient over 9 Gyr, 
and similar procedures for both the color gradients estimates 
and the relative analysis should be considered. In a forth- 
coming paper, taking into account all these factors, we will 
try to address the origin of the color gradients following their 
evolution from z^2 to z=0. 

For the remaining five ETGs of our sample, a pure radial 
variation of a single stellar population parameter is not able 
to reproduce the observed color gradients and global SEDs. 
Differently form the previous cases, where color gradients 
could be reproduced by a pure age or metallicity gradient 
as well as by a simultaneous radial variation of more than 
one parameter, these galaxies need a more complex scenario 
whereby more than one property of the stellar populations 
have to vary from the center to the periphery to generate 
the observed color gradients. 

Thus, our analysis clearly indicates that the properties 
of the stellar population and their distribution within high-z 
ETGs do not follow an homogeneous and common scheme. 

In the following we try to investigate if the theoreti- 
cal expectations of the widely accepted scenarios of galaxy 
formation can explain the observed color distributions. 

In the monolithic revised scenario, color gradients are 
supposed to be mainly due to metallicity gradient, being the 
contribution of age null or mild (and positive). Our findings 
of none correlation between color/metallicity gradients and 
total mass suggest that the monolithic revised scenario is 
not the favored mechanisms with which ETGs assembled 
their mass, although the narrow mass range covered and 
the assumption of the stellar mass as a proxy of the total 
stellar mass can affect this conclusion. 

Theoretical predictions of the inside- out- growth sce- 
nario point in f avor of compact ETGs with cores 
jWuvts et alJ koiO) redder than the outskirt regions. This 
negative color gradient seems to be due to a combined effect 
of negative metallicity a nd positive age gra dients, with a non 
negligible effect of dust (|Wuvts et al]|201Ch . To compare our 
results with the theoretical prediction of this model, we de- 
fine compact galaxies those ETGs with effective radius more 
than lcr smaller than those predicted by the local size-mass 
relation (SMR) for that mass. In Table 3 we report the com- 
pactness C for the galaxies of our sample defined as the ratio 
R e ,z=o/Re where R e is the effective radius of the galaxy and 
R e ,z=o is the radius that a galaxy of eq ual stellar mass w ould 
have at z=0 as derived by the SMR of lShen et ail (|2003l ). In 
our sample, compact galaxies, as we defined them, turn out 
to have C>2. In Fig [12] we report the SMR in the F850LP 
band for o ur sample (solid symbols) and for local galaxies 
(solid line, Shen et al. 2 0031 ). The dashed lines represent the 
scatter at la. It turns out that seven out of 11 galaxies (solid 



Spatially Resolved Colors and Stellar Population Properties in Early- Type Galaxies at z ~ 1.5. 



17 




Age [Gyr] Age [Gyr] 

Figure 11. Left panel: Black line shows the difference observed in the F850LP-F160W color of two stellar populations with age 
differentiating of 2Gyr, as a function of the age of the youngest stellar population. Red line shows the variation in the F850LP-F160W 
color observed in two populations with metallicities 2Zq and 0.2Zq, while blue line for two populations with metallicities Zq and 0.2Zq. 
Right panel Color gradients relative to the color variation in the left panel in the hypothesis that they occur between O.IRe and 3R e - 



points) of our sample are compact (magenta points) and 4 
are normal ETGs (cyan triangles). In Table 3, we report the 
classification in normal (N) and compact (C) for our galax- 
ies. Compact ETGs in our sample (as well as normal ETGs) 
show redder cores supporting the inside- out- growth scenario 
even if we cannot firmly establish the origin and the nature 
of this gradient. Actually, we cannot test the presence of a 
positive age gradient, since due to age/metallicity degener- 
acy we do not treat the case of simultaneous radial variation 
of age and metallicity. On the other hand, our results show 
that only two out of seven compact ETGs have color gradi- 
ents well reproduced by a metallicity gradient. 
Nonetheless, the simulations predicts that the negative color 
gradients produced by the interplay of age, metallicity and 
dust should result to cor relate with the integr ated rest-frame 
optical color. Following IWuvts et alJ |2010) in Fig. 1131 we 
show the ratio between the F850LP-band effective radius 
and F160W-band effective radius as a proxy of the color 
gradients versus the F850LP-F160W color. The points fol- 
low the same convention of Fig. 1121 We do not find any 
correlation neither in the whole sample, nor in the compact 
selection. The absence of such correlations cast some doubts 
on the validity of this scenario as a viable process to assem- 
ble the stellar mass in compact ETGs. 

Despite the fact that both the widely accepted formation 
scenarios do not seem to be able to reproduce the stel- 
lar content of high-z ETGs, we have to test the hypothesis 
whether ETGs can be assembled through a common forma- 
tion process and the distribution we observed in color gra- 
dients can be the final product of subsequent merger events. 
We have to bear in mind that our galaxies are all younger 
than 3.5 Gyr. This severely constrains the time each galaxy 
has at disposal to experience a merger event. Table 1 of 
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Figure 12. SM relation for local ETGs (solid line, Shen et al. 
2003) and for our sample (solid symbols). The dashed lines are 
the scatter lines at la. We shifted the Shen et al.s relation by 
a factor 1.2 towards lower masses to take into account the sys- 
tematic shift observed in the mass estimations using our models 
or those adopted by Shen et al. (2003). The circles are compact 
galaxies, that is, galaxies having the effective radius more than 
la smaller than those predicted by a local relation for that mass. 
On the contrary, galaxies having the effective radius comparable 
at la with those expected by Shen et al.s relations are classied as 
normal (triangle symbols). 
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Figure 13. The r e i?g50Lp/r e ,Fl60W ratio as proxy of F850LP- 
F160W color gradient vs the F850LP-F160W global color. 



Boyl an-Kolchin et all (|200cf ) shows the dynamical friction 
merging time in Gyr, r mer9er , for a host halo with virial 
mass Mh ost = 1O 12 M0, measured from numerical simula- 
tions. They assumed different initial orbital angular momen- 
tum, initial orbital energy, and ratio of initial satellite mass 
to initial host halo mass and only for two cases they found 
3.5 Gyr is enough to complete the merger, while in all the 
other cases T merger is greater than 4.3 Gyr. Thus, the differ- 
ent stellar content of high-z ETGs does not seem to be due 
to the effect of subsequent merger events, but primarily to 
the formation process. Actually, the continuum distribution 
of ETGs in the size- mass plane, both at high redshift and in 
the local Universe, together with the systematic direction of 
the color gradient (negative or null) of high-z ETGs, point 
toward a common formation process responsible of this con- 
tinuity. The possible different initial conditions, such as the 
different time scale of collapsing gas cloud, can be responsi- 
ble of the observed struc tural and dynamica l diffe rences as 
we previously suggested (|Saracco et al.ll201ll . |2012| V 
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Figure 7. (continued) 



Spatially Resolved Colors and Stellar Population Properties in Early- Type Galaxies at z ~ 1.5. 
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Figure 7. (continued) 



Spatially Resolved Colors and Stellar Population Properties in Early- Type Galaxies at z ~ 1.5. 
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Figure 10. Analysis of the radial variation of IMF's slope as possible driver of observed color gradients. Top and middle panels show the F850LP-F160W and F606W-F850LP color 
profiles, respectively (black lines), whit the relative fit (red lines). The bottom panel shows the observed global SED (black dots) as well as the stellar population parameters (black 
text) derived from its fit (black line) assuming a Salpeter IMF. Among all the stellar populations with varying IMF's slope but age, star-formation time scale, metallicity and dust 
extinction equal to the value obtained from the fit of the global SED, the red text identifies the one that best simultaneously reproduces the colors observed in the internal region (red 
dots in the top and middle panel). Similarly, blue text reports the value of the slope of the IMF that best reproduce the external colors (blue dots in the upper panels). The pure linear 
variation of the IMF's slope from the internal value to the external one, will produce the F850LP-F160W and F606W-F850LP color gradients represented with the dashed yellow lines. 
The contribution to the total stellar mass of the two populations is fixed in order to best fit the whole SED and are reported in the bottom panel (red and blue lines). 
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Figure 10. (continued) 
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